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Abstract 

The calculation of the SIR distribution at the typical receiver (or, equivalently, the success probability of 
transmissions over the typical link) in Poisson bipolar and cellular networks with Rayleigh fading is relatively 
straightforward, but it only provides limited information on the success probabilities of the individual links. 

This paper introduces the notion of the meta distribution of the SIR, which is the distribution of the conditional 
success probability P s given the point process, and provides bounds, an exact analytical expression, and a simple 
approximation for it. The meta distribution provides fine-grained information on the SIR and answers questions such 
as “What fraction of users in a Poisson cellular network achieve 90% link reliability if the required SIR is 5 dB?”. 

Interestingly, in the bipolar model, if the transmit probability p is reduced while increasing the network density 
A such that the density of concurrent transmitters A p stays constant as p —> 0, P s degenerates to a constant, i.e., 
all links have exactly the same success probability in the limit, which is the one of the typical link. In contrast, 
in the cellular case, if the interfering base stations are active independently with probability p, the variance of P s 
approaches a non-zero constant when p is reduced to 0 while keeping the mean success probability constant. 

Index Terms 

Stochastic geometry, Poisson point process, interference, SIR, coverage, cellular network, HetNets. 

I. Introduction 

A. Motivation 

Stochastic geometry provides the tools to analyze wireless networks with randomly placed nodes. A key quantity 
of interest in interference-limited networks is the success probability p s (6) = P(SIR > 6) of the transmission over 
the typical link, which corresponds to the complementary cumulative distribution (ccdf) of the signal-to-interference 
ratio (SIR). The calculation of p s involves spatial averaging, i.e., the evaluation of a certain expectation over the 
point process. While this expected value is certainly important, it does not reveal how concentrated the link success 
probabilities are. For example, in one network model, all links (or users) could have success probabilities between 
0.85 and 0.95, while in another, some links may have 0.5 and some may have 0.99. In both cases, we may find 
p s = 0.9, but the performances of the two networks in terms of connectivity, end-to-end delay, or quality-of- 
experience would differ greatly. Hence it is important to quantify the variability of the link reliabilities around 
Ps- 

To this end, our focus in this paper are random variables of the form 

P s {9) = P(SIR > 6 | <f>,tx), (1) 

where the conditional probability is taken over the fading and the channel access scheme (if random) of the 
interferes given the point process and given that the desired transmitter is active. The goal is to find (or bound) 
the ccdf of P s , defined as 

F Pe (x)±¥>'\P s (9)>x), x € [0,1], 


( 2 ) 
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where P !t denotes the reduced Palm measure of the point process, given that there is an active transmitter at the 
prescribed location. Since Fp a is the (complementary) distribution of a conditional probability, we call it the meta 
distribution of the SIR. Using this notation, the standard success probability is the mean 

p s (6) = E'-\P s (e)) = C FpSx) dx. 

Jo 

While a direct calculation of the ccdf (2) seems infeasible, we shall see that the moments of P s (9) can be expressed 
in closed-form, which allows the derivation of an exact analytical expression and simple bounds. The 6-th moment 
of P s {9) is denoted by Mb, i.e., we define 

M b {6) 4 E h (P s (9) b ) = C bx b ~ l F Pe {x) dx. 

Jo 

Hence we have p s (9) = M\{9). 

B. Contributions 

The contributions of the paper are: 

• We introduce the meta distribution of the SIR. 

• We give closed-form expression of the moments M b for Poisson bipolar networks with ALOHA and for 
Poisson cellular networks, both for Rayleigh fading. 

• We provide an analytical expression for the exact meta distribution for the two types of networks. 

• We propose the beta distribution as a highly accurate approximation. 

• We show that, remarkably, in the limit of very dense bipolar networks with small transmit probability, all links 
have the same success probability. This is not the case in cellular networks with random (interfering) base 
station activity, since the variance M 2 — Mf is bounded away from zero when the probability of a base station 
being active goes to 0. 

• We give the conditions on the SIR threshold 9 and the transmit probability p for a finite mean local delay. 

C. Related work 

The calculation of the (mean) success probability p s (9) in Poisson bipolar networks is provided in [1] but can 
be traced back to [2], In [3], the moments Mb of the link success probabilities are calculated under the assumption 
of no MAC scheme (i.e., all nodes always transmit), and bounds on the distribution are obtained. 

For Poisson cellular models, where the typical user is associated with the nearest base station (strongest base 
station on average), the result was derived in [4] and extended to the multi-tier Poisson case (HIP model) in [5]. 

The joint success probability of multiple transmissions in Poisson bipolar networks is calculated in [6]. Similarly, 
[7] determined the joint success probabilities of multiple transmissions (or transmissions over multiple resource 
blocks) for Poisson cellular networks. As we shall see, these joint probabilities are related to the integer moments 
M;,, of the conditional success probabilities. 

D. The meta distribution 

In this section, we formally introduce the concept of a meta distribution, which is the distribution of the conditional 
distribution P s . 

Definition 1 (Meta distribution) The meta distribution of the SIR is the two-parameter distribution function 

F(9,x) = F Pe (9, x) = P !t (P s (0) >x), 9e M + , x € [0,1]. 
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We have F(0,x) = 1 for x < 1, lim^oo F(9, x) = 0 for x > 0, F(Q, 0) = 1, and F(6, 1) = 0. For fixed 9, it 
is a standard ccdf and yields the probability that the typical link or user achieves an SIR of 9 or, equivalently, 
the fraction of links or users (assuming a uniform user distribution) that achieve this SIR. Generally, it yields the 
fraction of links or users that achieve an SIR of 9 with probability at least x. 

In the next two sections, we will calculate the meta distribution and bounds for Poisson bipolar and cellular 
networks, respectively. 


II. Poisson Bipolar Networks 

A. System Model 

We consider the Poisson bipolar model [8, Def. 5.8], where the (potential) transmitters form a Poisson point 
process (PPP) ( I> of intensity A and each one has a dedicated receiver at distance R in a random orientation. In each 
time slot, nodes in <3? independently transmit with probability p, and all channels are subject to Rayleigh fading. 

We use the standard path loss model with exponent a, define 6 = 2/a, and we let C = A7rf? 2 r(l — b')r(l + 8) 
be a coefficient that does not depend on 9. The success probability of the typical link is well known, see, e.g., [1], 
[8], [9], and can be expressed as 

Ps (9) 4 P !t (SIR > 9) = Afi(0) = e~ ce6p . 

Due to the ergodicity of the PPP, the ccdf of P s can be alternatively written as the limit 

F Pe (x) = lim 1 l(P(SIRy > 9 | $) > x), 

i—roo A»7rr z z ' 
y£<t> 

\\y\\< r 

where y is the receiver of transmitter y and l(-) is the indicator function. This shows that Fp s (x) denotes the 
fraction of links in the network (in each realization of the point process) that, when scheduled to transmit 1 , have 
a success probability larger than x. 

The link success probabilities for a given realization can also be “attached” to each point of the transmitter 
process T to form a marked point process <I> = {(jp. P/')}- The meta distribution can then be interpreted as 
the mark distribution, parametrized by 9. Due to the interference correlation [10], the marks of nearby nodes are 
correlated, hence is not an independently marked process. 

Fig. 1 shows an example realization of a Poisson bipolar network together with the success probabilities for 
each link, averaged over the fading and ALOHA. As expected, links whose receivers are relatively isolated from 
interfering transmitters have a high success rate, while those in crowded parts of the network suffer from a low 
one. 


B. Moments 
Let 

For p = 1, 


D b (p,S)±J2 


k =1 


b\f8-l 


k \k- 1 


D b (l,S) = 


p k , be C and p, 5 € [0,1]. 


r(6 + <y) 


(3) 


r(&)r(i + <$)’ 

which is not defined if b G IF or b + 8 G Z - . For 8 € {0,1}, the function simplifies to D b (p, 0) = 1 — (1 — p) b 
and D b (p, 1) = bp. 

'The received signal power is assumed zero if the desired transmitter is not active, so the SIR is zero in this case. 
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Fig. 1. Realization of a Poisson bipolar network for A = 1, R = 1/2, p = 1/2, 9 = 1, a = 4, resulting in p s = 0.54. The number next to 
each link is its success probability (averaged over fading and ALOHA). 


Alternatively, the function can be expressed using the Gaussian hypergeometric function y P\ as 

Db(p,S) =pb 2 F 1 (l - b, 1 - 8;2;p). (4) 

Theorem 1 (Moments for bipolar network with ALOHA) Given that the typical link is active, the moment Mb 
of the conditional success probability is 

M b [9) = exp (-Ce s D b (p, <5)) , b€ C, (5) 

whenever Db(p , 6) is defined. 

Proof: See Appendix A. 

An important and helpful observation in the proof is that the calculation of the n-th moment for n € N is the 
same as that of the joint success probability of n transmissions, calculated in [6]. In this case, D n (p,6) is given 
by the finite sum 

n 

D n (p,5) = ^2 

k =1 

which is a polynomial of degree n in p and degree n — 1 in <5 and called the diversity polynomial in [6, Def. 1], 
Since (5) is valid for (essentially) any b G C, we can use it to obtain the -1-st moment as 

M_i(0) = exp(CO s p(l — p) 5-1 ) 

= M“ (1 “ p) *"\ p<l. (6) 

M- 1 is the mean number of transmission attempts needed to succeed once if the transmitter is allowed to keep 
transmitting until success. This quantity is termed mean local delay and is calculated in [11, Lemma 2]. Noteworthy 
is the phase transition at p = 1. For p = 1 — e, the mean local delay is finite for all e > 0. But if all nodes always 
transmit, it is infinite. 


<5—1 
k - 1 


P 
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An interesting question is what happens when p —> 0 while the transmitter density p\ (and thus M\) is kept 
constant. It is answered in the following corollary. 


Corollary 1 (Concentration as p —> 0) Denoting the transmitter density as r = A p and keeping it (and thus Mi) 
fixed while letting p -A 0, we have 

lim P s {6) = p s {0) 

p—>■ o 

A P=T 

in mean square (and probability and distribution). 

Proof: From (5), the second moment is 

M 2 (6) = e -^(2 P +(«-i) P a ) ) 


and the variance, expressed in terms of M\ (which is kept constant), is 

var P s (9) = - 1). 


It follows that 


lim vax PJO) = 0. 

p —^0 

A p=r 


(7) 


So if CO s p is kept constant, the vaiiance can be adjusted by changing p. For example, if C = 1/(10 pd s ), 
Mi = e -1 / 10 ~ 0.9, and the variance can be reduced to 0 by letting p -A 0. So, counterintuitively, a small p 
decreases the vaiiance and, in the limit, all links in the network have exactly the same success probability. 

More precisely, the variance is proportional to p for small p if M\ is kept constant: 

v&t P s (9) ~ —Mi log(Mi)(l — 5)p, p —> 0. 

The next result provides tight bounds on the moments if p = 1 for b € M + . and indicate upper bound 


and lower bounds with asymptotic equality (here as 6 —>• oo), respectively. 


Corollary 2 (Bounds on moments for p = 1) For b > 0, 

r(b+«) 

M b = AT 1 r(1+5)r(6) > exp(—C9 s b s ), 

(8) 

for b > 1, 

M b < M b i, 

(9) 

and for 0 < 6 < 1, 

M b > Mf . 

(10) 


Proof: The lower bound (8) follows from (5) by setting p = 1 and the asymptotic bound T(b + 5)/T(b) < if 
for b > 0. Conversely, T(7; + A) /T(b) > If Y (1 + A) for all b > 1, which yields the upper bound (9): 

M h < exp(—C6 5 T(1 + A)) = Aff, b > 1. 

For b < 1, T(b + A) /T(b) < lfT(l + A), and the direction of the inequality is reversed, yielding (10). ■ 

The third bound is tighter than the first one in the regime where it is valid. Further, since 

Mf = exp C(b6 , 

the 6-th moment is bounded by the first moment evaluated at b6, i.e., 

M b (0) < Mi(b0), 6 > 1 , 


and vice versa if b < 1. 
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C. Exact expression 

An exact integral expression can be obtained from the purely imaginary moments Mjt, t € R, j = \J— 1. 
Corollary 3 (Exact integral expression) The meta distribution is given by 

1 ncr nr -U noSpsmx, 

-df, 


. 1 1 e CdS ^i D it) sin(f logx + C9 s ^s(Dp)) 

F(e ’ X)= 2 -*J 0 - t - 


( 11 ) 


where Dj t = Dj t (p,5) is given in (3) and 9?(z) and 9?(z) denote the real and imaginary parts of the complex 
number z, respectively. 


Proof: Let X = log P s {9). The characteristic function of X i 


is 


<px(t) = Ee jtx = E (P s (6f) = M jt , t G R. 

where M ]t is given in (5). Then by the Gil-Pelaez theorem [12], the ccdf of X is given by 

D Q(e~i tx M jt ) 


- , . 1 1 
F x(x) = - + - 

Z 7r 


t 


-df. 


( 12 ) 


(13) 


Since P (P s {0) > x) = P(log.P s (#) > logx), 

— . , 1 1 f°° cj( e -f* lo g x M jt ) , 

Fp. (x) = - + ~ - ---—df, 

2 ir J 0 t. 

and the result follows from Thm. 1 and some simplification. ■ 

Since \M jt \ essentially decreases exponentially with f, this integral can be evaluated very efficiently. The curve 
marked with o in Fig. 2 shows the exact meta distribution F(l,x) for A p = 1/4 with different values of A and 
p. As predicted by Cor. 1, the variance of P s is reduced when p is smaller. Next we will derive the bounds also 
shown in the figure. 


D. Classical bounds on the meta distribution 

Simple bounds on the meta distribution can be established using classical methods. 

Corollary 4 (Markov and Chebyshev bounds) For x G [0,1], the meta distribution is bounded as 


(i — x y 


Let V = var P s {0) = M -2 — Mf. For x < M±, 


while for x > M\, 

Lastly, 


Fp e (x) > 1 - 
FpAx) < 


V 


(x — M\) 2 ’ 
V 


F Ps (xMi) > 


(x — M\) 2 ' 
(1 — x) 2 


ie (o,i). 


(14) 


(15) 


(16) 


(17) 


1-Mf {1 ^^(l-x) 2 ’ 

Proof: (14) follows from Markov’s inequality, while (15) and (16) follow from Chebyshev’s inequality. The 
lower bound (17) is the Paley-Zygmund (or Cauchy-Schwarz) bound. ■ 

For the lower (or reverse) Markov bound in (14), the integer moments of 1 — P s {0) are easily found using binomial 
expansion. For b = — 1, the Markov inequality also yields the lower bound Fp s (x) > 1 — xM_i, where M_i is 
given in (6). 
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(a) A = 1 ,p = 1/4, and var(P s ) = 0.0212. (b) A = 5, p = 1/20, and var(P s ) = 0.00418 

Fig. 2. The exact meta distribution (11) and the Markov bounds (14) for b £ [4], (15), and (16) for a = 4, 6 = 1. R = 1/2, and A p = 1/4. 
The resulting mean success probability is p 3 = Mi = 0.735. The variance depends on the values of p and A; it is proportional to p for 
small p. 


These bounds are illustrated in the two plots in Fig. 2. For the Markov bounds, the four lower and upper bounds 
correspond to b = 1,2,3,4. It is apparent that the variance decreases with decreasing p and that the bounds get 
tighter also. 

Written differently, (15) and (16) state that 

M s ~ 1 — 1 

Fp s (q M i) >1 77-3^n, 0 < q < 1, 


(1 - 4) 5 


and 


M 5-1 — 1 

Fp B {qM\) < ~7T - ZyT- L < T < M-\ 


(1 ~q) 2 

The upper bound is useful for small Mi, while the lower bound is useful for Mi ss 1. 

So as p —» 0, F(P s (0) > xM\) —>• 1 Vx € (0,1), in accordance with Cor. 1. 

The Paley-Zygmund bound is useful to bound the fraction of links that has at least a certain fraction of the 
average performance. For example, the fraction of links having better than half the average reliability is lower 
bounded as 


F h (P s (9) > Mi/2) > 


1/4 


5/4-Mf (1 ^ 

As p —>• 0, the lower bound approaches 1, again as expected from the concentration result in Cor. 1. 


E. Best bounds given four moments 

Here we establish the tightest possible lower and upper bounds on the distribution given the first four moments. 
Generally, this problem can be formulated as follows. Letting be the class of distributions (cdfs) with moments 
Mi,...,Mfc, we would like to find 

L(x) = min F(x), x € (0,1) 


and 


F£M k 


U(x) = max F(x), x € (0,1). 
FeMk 

















So for each x in the support of the distribution, we would like to find the minimum and maximum over all 
distributions with the prescribed k moments. To find L and U for k = 4, we are applying the method from [13]. 
It determines the best lower and upper bounds 


L(x) < F y (x) < U(x) 


given the four moments E(Y k ), k € [4], for a general continuous random variable Y. 

To bound the cdf Fy(x) at a target value x, first the moments are calculated for the random variable shifted by 
x so that the new target location is 0, i.e., 

rrii(x) = f (y - x ) 1 dF Y (y) 

J o 


J2[l)(-xr k HY k ), x € [0,1]. 


Using these shifted means, following [13], we define (omitting the dependence on x of the shifted moments to 
avoid overly cumbrous notation) 

q{x ) = \J (—m 2 rri 3 + — 4 (m| — — 7772777 , 4 ) 

a ~ m 2 + 2 mim. 2 TO 3 — — 777(7774 + 77727774 


Po(x) = 


7772777,4 — 


y\{x) = 


A 77127713 — 77117714 — q{x) 


2(m| — m\77l3) 


V2{x) = 


A 77127713 — 77117714 + q{x) 


2(7712 ~ mi77l3) 


A rn%- 77747713 f 
n(x) = — (■’ 
Pl{x) = 1 ~Po(x) ~P 2 (x), 


('7772 — 277117712777,3 + TTl\m 4) (—77727773 + 77747774 + q(x) 

2(7772 ~ 77l47773)(—77l^ + 77727774) 


and the bounds follow as 


1 Pi(x)+p2(x) if yi(x) < 0, y 2 (x) < 0 
Pi(x) if yi(x) < 0, 2 / 2 0) > 0 

0 if yi(x) > 0 , y 2 (x) > 0 

1 if 2/i (x) < 0, r/2(x) < 0 

U(x) = < Po (x) +p4(x) if 2 /i (x) < 0 , 2/2(x) > 0 

,Po(x) if 2/1 (x) > o, 2/2(x) > 0 

Since q(x) > 0, it is not possible that 2/1 (x) > 0 and 2/2(x) < 0. 

In our application Y = P s (9 ), E(Y k ) = M&, and since we are working with ccdfs, we have 


1 — U(x) < F(6,x) < 1 — L(x). 


Fig. 3 shows these best bounds, together with the lower and upper envelopes of the Markov upper and lower 
bounds for b € [4] and the Paley-Zygmund lower bound. In some intervals, the classical bounds are near-optimum, 
while in others, the best bounds are significantly tighter. 

The method in [13] is not restricted to four moments, but it is considerably more tedious to apply if more 
moments are considered. 





(a) A = 1 => p s = 0.54, var(P s ) = 0.049. (b) A = 1/5 => p s = 0.88, var(P s ) = 0.024 

Fig. 3. The exact meta distribution (11), the best Markov bounds (14) for b £ [4], and the best overall bounds per (18) and (19) (given the 
first four moments) for a = 4, 9 = 1, R = 1/2, and p = 1/2. The reduction of A from 1 to 1/5 results in a reduction of the variance of 
only 1/2, since p stays the same. 



k = -1 

CO 

II 

k = 4 

k = 5 

k — 6 

k = 7 

OO 

II 

M k 

1.4278 

0.4418 

0.3571 

0.2947 

0.2476 

0.2110 

0.1820 

E(X k ) 

1.4333 

0.4412 

0.3555 

0.2921 

0.2440 

0.2066 

0.1770 

ratio 

0.9962 

1.0014 

1.0044 

1.0090 

1.0147 

1.0211 

1.0280 


TABLE I 

Comparison of moments Mb and E(X k ) of the beta approximation for the parameter set in Fig. 2(a). 


F. Approximation with beta distribution 


Since P s (0) is supported on [0,1], a natural choice for a simple approximating distribution is the beta distribution. 
The probability density function (pdf) of a beta distributed random variable X with mean p is 

- , V _ X (1 - x)P 1 

fx[X) - B(p(3 /(1 — p), /3) ’ 

where B(-, •) is the beta function. The variance is given by 


2 A y 

a = var X 


m(1 - f) 2 
(3 + 1 — p 


Matching mean and variance a 2 yields p = Mi and 

(3 _/r(l-/r) 2 n ^ _ (p - M 2 )(l - p) 

P a 2 1 N M 2 ~p 2 ' 

As illustrated in Fig. 4 (same parameters as in Figs. 2 and 3), the beta distribution provides an excellent match 
for the distribution of the link success probabilities, which is also corroborated by the fact that the higher moments 
E(X fc ) of the matched beta distribution are very close to M*.. For example, for the parameters in Fig. 2(a), the 
analytical -1-st and 3-rd through 8-th moments differ by less than 3%, as shown in Table I. So the skewness and 
kurtosis and the mean local delay are approximated very accurately also. 




































10 




0.9 

0.8 

0.7 

0.6 

0.5 

0.4 

0.3 

0.2 

0.1 

0 


-e- exact 

— beta approximation 
0.2 0.4 


VltmtWnnim EBBB8886BBB8aaBMB 


4 

0.9 

0.8 

0.7 

0.6 

0.5 

0.4 

0.3 

0.2 

0.1 

°0 


-e- exact 

— beta approximation 
0.2 0.4 


(a) Parameters from Fig. 2 (a) and (b). (b) Parameters from Fig. 3 (a) and (b). 

Fig. 4. The exact meta distribution and the beta distribution approximation for the two sets of parameters considered in the plots of Figs. 2 
and 3. 



Fig. 5. Three-dimensional plot of the meta distribution F(9,x) for A = 1, p — 1/4, a = 4, and R = 1/2. 


G. Illustrations of the meta distribution 

An illustration of the meta distribution is shown in Fig. 5. It shows qualitatively that, for the chosen parameters, 
most links achieve an SIR of —10 dB with probability 80%, while an SIR of 10 is achieved with probability 80% 
by virtually no links. For quantitative purposes, the cross-sections and contours are more informative, as shown in 
the next figures. 

Fig. 6(a) enables a more precise statement about the fraction of links achieving an SIR of —10 dB with 80% 
reliability—it is 0.93. It also shows that at 6 = 0 dB, 60% of the links have a success probability of at least 80%. 

As a function of 6 for fixed x, the value of 6 can be determined such that at least a fraction x of users have a 
success probability p m ; n . For example, Fig. 6(b) shows that to achieve at least 80% success probability for 80% of 
the links, a 6 of at most —7.6 dB can be chosen. 

The contour plot Fig. 7 visualizes the trade-off between x and 9. It shows the combinations (9, x) that can be 
achieved by a certain fraction of links u. For example, the curve for link fraction u = 0.95 shows that 95% of the 
links achieve an SIR of —5 dB with probability 0.6 and an SIR of 5 dB with probability 0.31. 
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(a) Meta distribution for 9 = -10,-5,0,5,10,15 dB. The curve (b) Meta distribution as a function of 9 for 
for 9 = 0 dB is marked with o. x = 0.4, 0.5, 0.6, 0.7, 0.8, 0.9. 

Fig. 6. Cross-sections through the meta distribution along the x and 9 axes for A = 1, p = 1/4, a = 4, R = 1/2. 



Fig. 7. Contour plot of meta distribution F(9, x) for A = 1, p = 1/4, a = 4, and R = 1/2. The values at the curves are F(6, x) = u = 
0.5, 0.6, 0.7, 0.8, 0.9, 0.95 (from top to bottom). 


Hence the contour plot illustrates and quantities the trade-off between data rate (as determined by 0) and reliability 
(given by the parameter x ) in bipolar networks. 

III. Poisson Cellular Networks 

A. System model 

In Poisson cellular networks, base stations (BSs) form a PPP of intensity A, while users form a stationary point 
process of intensity A u . We focus on the downlink and on nearest-BS association, /.<?., each BS serves all the users 
in its Voronoi cell, and first assume that all BSs are always active. An example realization where users form a 
square lattice is shown in Fig. 8. 
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Fig. 8. Realization of a Poisson cellular network with BS density A = 1, users forming a square lattice of density A u = 3, d = 1, and 
a = 4, resulting in p B = 0.56. The BSs are indicated by x and the users by o. The number next to each user is its success probability 
(averaged over fading) or its mark, and the dashed lines are the edges of the Voronoi cells of the BS PPP. 


As in the bipolar case, we assume the standard path loss law with path loss exponent a = 2/5 and Rayleigh 
fading. The standard (mean) success probability (or SIR distribution) is the success probability of the typical user, 
assumed at the origin o, which is known from [4] as 


Ps(0) = P°(SIR > 9) = 


1 


2 *i(1, - 6 ] l- <5; - 0 )' 

The probability also has a spatial interpretation: for each realization of the BS and user point processes, it gives 
the fraction of users achieving an SIR of at least 9 in a given time slot. It depends neither on the user density nor 
on the BS density. 

Again we define the conditional success probability 


P s {9) = P°(SIR >9 | $), 


which is the probability that the SIR at the origin exceeds 9 given the BS process and given that a user is located 
at o. The quantity of interest is the meta distribution of the SIR, which is the distribution (ccdf) of P s : 

F(9,x) = Fp s (x) =P(P s (0) >x), 9£R + , x€ [0,1] 

It gives detailed information about the user experience by providing the fraction of users achieving an SIR of 9 
with reliability at least x. 

As before, a direct calculation of this meta distribution seems infeasible and we thus focus on the moments 

M b = E(P s (9) b ) first. 


B. Moments 


Theorem 2 (Moments for cellular network) The moments of the conditional success probability for Poisson 
cellular networks are given by 


M b = 


1 


2 Fi(b,-5-i-5--ey 


b e 


( 20 ) 
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Proof: Let xo = arg min{x € 4>: ||x||} be the serving BS of the typical user. Given the BS process <I\ the 
success probability is 


P,(6) = v(h > 


n 


o:G$\{a: 0 } 


*0 22 h x \\ x \ 

a;G$\{a;o} 

1 

l + 0(||xo||/||x||)“- 


4* 


The 6-th moment follows as 


M b = E ] [ 




(1 + 0 ( 11 x 011 / 11 x 11 )^- 


( 21 ) 


Instead of calculating this expectation in two steps as usual (first condition on ||xo|| then take the expectation 
w.r.t. it), we use the recent result [14, Lemma 1], which requires the calculation of only one finite integral. The 
lemma gives the pgfl of the relative distance process (RDP ), defined as 


77 = {x £ 4> \ {x 0 } : ||x 0 ||/||x||}, 


when 4> is a PPP. Since (21), depends on the BS locations only through the relative distances, we can directly 
apply the pgfl of the RDP and obtain 


M b 


1 + 2/ (l 

o v 


1 

(l+Qr+O r_3dr 


( 22 ) 


which can be expressed as (20). ■ 

Sometimes the calculation of the hypergeometric function with negative last argument can cause numerical 
problems. In such cases, the alternative form 


M b = 


(1 + 0) b 


2 Fi( 6,1;1-6;0/(l+ 0))’ 

obtained through Euler’s transformation, is helpful. 

For 6 = — 1, (20) (or (22)—no “detour” using hypergeometric functions needed in this case) simplifies to 

1-6 


M_ i = 


1 - 6(1 + 0 )’ 


0 < 1 / 6 - 1 . 


(23) 


As in the bipolar case, this is the mean local delay if 0 < 1/6 — 1. Converseley, if 0 > a/2 — 1, the mean local 
delay is infinite due to the correlated interference in the system. This phase transition in the mean local delay is 
similar to the one observed in [6], [11], [15] for ad hoc networks. Incidentally, the condition can also be expressed 
as 0 MISR < 1, where MISR is the mean interference-to-signal ratio of the PPP introduced in [16]. 

For 6 € N, the moment M b equals the joint success probability of 6 transmissions, which was calculated in [7, 
Thm. 2] using a different (less direct) method. 

Fig. 9 shows the standard success probability M\ = p s and the variance as a function of 0 for a = 3,4. Since 
the variance necessarily tends to zero for both 0->O and 0 —>• oo, it assumes a maximum at some finite value of 
0. A numerical evaluation shows that for a = 3, the variance is maximized quite exactly at 0 = 1, and for both 
values of a, the success probability at which the variance is maximized is p s = 0.38. 
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Fig. 9. Success probability Mi and variance M 2 — Alf for a = 3 and a = 4. 




(a) 9 = 1 => p a = 0.56, var(Ps) = 0.098 (b) 9 = 1/10 =4> p B = 0.91, var(P s ) = 0.0086 

Fig. 10. The exact meta distribution (24), the best Markov bounds (14) for b £ [4], the Paley-Zygmund lower bound, and the best overall 
bounds (given the first four moments) for a = 4. 


C. Exact expression, bounds, and beta approximation 

As in the bipolar case, we obtain an exact expression for the meta distribution from the Gil-Pelaez theorem. 


Corollary 5 The SIR meta distribution for Poisson cellular networks is given by 


F{6,x) = \ + - 

A 7T 




d t 


(24) 


Numerical investigations indicate that \Mj t \ = 0(£ -1 ), t —> 00, so the integrand decays with t~ 2 and the integral 
can be evaluated efficiently. 

Fig. 10 shows the exact distribution and the classical and best bounds for 6 = 1 and 6 = 1/10, respectively. 
Interestingly, the meta distribution F( 1, x) has almost constant slope, which means that the user success probabilities 
are essentially uniformly distributed between 0 and 1. 
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Fig. 11. Exact ccdf and beta approximation for 9 = 1/10,1,10 for a = 4. 



Fig. 12. Contour plot of meta distribution F(9, x) for a = 4. The values at the curves are F(9, x) = u = 0.5, 0.6, 0.7, 0.8, 0.9, 0.95 (from 
top to bottom). 


Fig. 11 shows that the beta approximation provides an excellent fit over a wide range of 6 values. It also serves 
as an illustration of the meta distribution showing what combinations of reliability x and fraction of users can be 
achieved for 9 € {—10,0,10} dB. 

Lastly, Fig. 12 shows a contour plot of the meta distribution for a = 4. An operator who is interested in the 
performance of the “5% user” (the user in the bottom 5-th percentile in terms of performance) can use the bottom 
curve, corresponding to F(6, x ) = 0.95, to find the performance trade-off that such a user can achieve. For example, 
it can achieve an SIR of —10 dB with reliability 0.72 or an SIR of —4.3 dB with reliability 0.3. 

D. Effect of random base station activity 

Flere we investigate the effect on the meta distribution if interfering BSs were active only with probability p. 
This is similar to the model studied in [4, Sec. VI], where a frequency reuse parameter n was introduced and each 
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BS is assumed to choose one of k bands independently at random. Hence the two models are the same if we set 
p = n 1 (apart from the fact that k 6 N, whereas no such restriction is imposed on p '). 

Theorem 3 The b-th moment of the success probability in a Poisson cellular network where interfering BSs are 
active independently with probability p can be expressed as 

-l 


M b (p) = ( 1 - jr (tj) (- P 0) k 2^1 (*, k-5'k + l-5--9) 

\ k =l ' ' / 


(25) 


1 + 9r a 


+ 1 —p 


Proof: If interfering BSs are active independently with probability p in each time slot, we have 

(~ P 

r£TZ 

and thus 

pOr° 


p s (o) = n 


M h (p )=e n i 


re TZ 


1 + 9r a 


Hence we need to modify (22) to 


M b (p) = 


(26) 


1 + 2/ 1-1- 


pOr 01 
1 +0 r ° 


3 dr 


For general b G C, letting x = r a , the integral in (26) can be expanded as 2 


E 

fc=i 


b\-(- P 9) k + 
k 


a 


q \ 1 + Ox 


x 5 1 d.x = 


E 

fc=i 


b\ ~(~p9) h 

k) ka — 2 


2 F 1 {k,k-6;k + l-6-,-0), (27) 


and we obtain the result. 

For b = 1, this yields the success probability 

Ps{0,p) = 


1 


1 + P@t^$ 2+i (1) 1 — <5; 2 — (5, -9) 

1 


l -p + P2+i(l, -5; l - <5; - 0 ) 

The first expression corresponds to [4, Eqn. (19)], while the second one follows from the identity 

95 


1-5 

For b = — 1, (26) yields 

M_! = 


2 +i(1 1 1 - <5; 2 - 5; -9) + 1 = 2 Fi(1, -<5; 1 - <5; -9). 


P < PciQ). 


(28) 

(29) 

(30) 

(31) 


1 - pOj^s 2^i(l. 1 - -5; 2 - 5, -9(1 - p))' 

Here p c (9) is the critical transmit probability denoting the phase transition from finite to infinite mean local delay. 
If 9 < 1/5 — 1, we know from (23) that p c {9) = 1. If p < 1, a larger 9 can be accommodated while maintaining 
a finite mean local delay. Fig. 13 shows the critical probability p c (9) and two conjectured bounds, which arc 

Pc(0) > (j^9)- s / 2 and p c {9) < ( T ^9)~ s . 

Next we provide an asymptotic result on the success probability p s (p, 9) as p -x 0 while keeping p9 s constant. 


2 See the appendix, where a similar technique is used. 
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Fig. 13. Critical probability p c (in dB) for finite mean local delay as a function of 6 for a = 3,4 and conjectured lower and upper bounds. 



Fig. 14. Contour plot showing the combinations of 8 and p (in dB) that achieve a given target success probability p t G 
{0.1, 0.2,..., 0.9, 0.95} for a = 4. The dashed lines are the asymptotes obtained from (32). 

Corollary 6 Let t = p9 s . As p ^ 0 and 0 —>■ oc such that t stays constant, 

p s {6,p) ~ 1 } 5/ . —t = ( 32 ) 

I+ptJ°/smcd c + sined 

Proof: From Thm. 4 and Lemma 6 in [14], 2 -F) (1, —<5; 1 — 5; —6 ) ~ 9 s / sine 5, 9 —>• oo. Inserting this in (29) 
and letting p —>• 0 and 0 oc while keeping p9 s constant yields the result. ■ 

The corollary implies that 

Ps(0,p) ~ P s (c 1/5 0,p/c), c > 1. 

So in the limit of small p, if p is decreased by 10 dB, 9 can be increased by 5 a dB to maintain the same success 
probability. 

Fig. 14 shows a contour plot indicating the combinations of 8 and p (in dB) that achieve a given target success 
probability p t , together with the asymptotes obtained from (32) by calculating t from t = (j\ 1 — 1) sine 6 and then 
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Fig. 15. Variance M 2 — Mi as a function of the BS activity probability p for target success probabilities p t £ {0.7, 0.8, 0.9} for a = 4. 
The dashed lines are the asymptotes from (34). 


plotting 9(jp) = (f/p) l /V \ which is a line in the log-log plot. Hence, keeping p6 5 constant results asymptotically in 
the same success probability, as p —>■ 0 or 6 —>• 00; in contrast, in the bipolar case, keeping p9 s constant results in 
exacty the same success probability for all values of p and 6. 

An important question is whether—as in the bipolar case—the variance goes to 0 as p —>• 0 while keeping p s 
constant. The last corollary answers that question. 


Corollary 7 Given t = p9°, 

lim varP s (6*,p) = '^ 11C ^ _ 

p— to y 2t + sine <5 

e=(t/ P y /' 5 

Expressed as a function of the target success probability p t , 

lim var P s (9,p) = ——— 

p-s>o 2 — p t 

e=(t/p) 1/s 


sine S 
t + sine 5 


Pi 


(33) 


(34) 


Proof: The inverse of the second moment follows from Thm. 3 and is given by 

Mf 1 = 1 + 2 p 9 -A- 2 Fi (1, 1-5-2-6,-0)- p 2 9 2 -^— 2 +(2, 2-6; 3-6,-0). 
i—o 2—o 

^ ---/ ^' 

A B 

As 9 -A- 00, combining (32) and (28), A = 9 s / sine 6. For B, we have 3 B = 0(0 <5 ). Hence, for some constant 
c > 0, 

lim Mf 1 = 1 + 2t/ sine <5 — ptc = 1 + 2 1/ sine 6. 
p—>- 0 

o=(t/p) 1/s 


The result follows from var+ s = M 2 — Mf, with M\ given in (32). ■ 

Fig. 15 displays the variance as a function of p for different target success probabilities. These are the variances 
obtained along the corresponding contour lines in Fig. 14. The asymptotic variance from (34) is also shown. It 
can be seen that the transmit probability has relatively little impact on the variance, especially for higher success 
probabilities. So, in contrast to the bipolar case, the disparity in the user experience cannot be significantly reduced 
by random BS activation patterns. 


3 See, e.g., http://dlmf.nist.gOv/15.8#E2. 
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IV. Conclusions 

While spatial averages, such as the success probability of a transmission over the typical link (or standard SIR 
distribution), arc useful, they do not provide much information about the performance of the individual links or 
users in a given realization of the network. To overcome this drawback, this paper introduces the meta distribution 
of the SIR, which is the distribution of the conditional SIR distribution (or success probability) given the point 
process, and provides an exact expression, bounds, and an approximation, for Poisson bipolar and cellular networks. 
Hence the complete distribution of the conditional link success probability P s in both types of Poisson networks 
can be characterized. The complete distribution of PJ9) provides much more fine-grained information that just the 
mean p s (9) that is usually consiered. 

The key insight is that the moments of P s can be calculated in closed-form. Hence standard and optimum 
moment-based bounding techniques can be employed, which yield lower and upper bounds that arc reasonably 
tight in some regimes. Moreover, an approximation by a beta distribution by matching first and second moments 
turns out to be matching the exact distributions extremely accurately. 

Bipolar networks with ALOHA exhibit the interesting property that the variance of P s goes to 0 as the transmit 
probability p —»• 0 while keeping the (mean) success probability constant. This is, however, not the case for cellular 
networks. If interfering base stations are active independently with probability p, the variance approaches a non-zero 
constant as p —>• 0, again while keeping a constant success probability p a . So the deployment of an ultra-dense 
network of small cells that are only active with small probability (when a user requires service in their cell) does 
not significantly reduce the disparity of user experiences. On the positive side, lowering p allows an increase of 9 
without affecting p s . To be precise, decreasing p by 10 dB allows an increase of 6 by 5a dB. 

From a broader perspective, the results show that it is possible in certain cases to not only derive spatial averages, 
but complete spatial distributions, which constitute rather shaip results on the network performance since they 
capture the statistics of all links in a given realization of the network. Hence it is demonstrated that stochastic 
geometry allows for the calculation of (even) stronger results than spatial averages. 
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Appendix 

A. Proof of Theorem 1 

Proof: Given <I>, the success probability is 

P s (9) = P (h > 9'I | $) = E(e- 0 ' 7 | $), 

where 6' = 9R a and 

I = h x \\x\\~ a l{x € 4>t). 




Averaging over the fading and ALOHA, it follows that 


«(«) = n y 


p 




+ 9'\\x\ 


+ 1 — p. 
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Hence we have 


M b = E 


n 


P 


Lxe$ 

= exp (—A 


1 + #'11x11 
1 - 


+ 1 —p 


P 

l + #'||x| 


1 — p 


dx 


This is the same integral as in [6, Appendix A] and thus for 6 £ N, the resulting expression is the diversity 
polynomial derived there. 

For general (non-integer) b, the proof in [6, Appendix A] needs to be modified. Expressing the moments as 
M b = e A Fb , we have from (29) in that paper 

3/ \ bl 


F b = i t5 


1 - 1 - 


p6 


u + 6' 


u s ~ l du. 


For general 6 G C, we replace the summation bound by oo since 


and we obtain 


a - = E J 


k =0 


r oo 00 


F b = ir6 £ (-1 ) k+ HpO') k - 


“i 


,5-1 


(u + 9') k 


du 


k=\ 


I o (u + 6'Y 


d u. 


For the integral we have 


and thus 


"°° u5 ~ l - du = d /5 - k{ - 1)k+ln r(5) 


/ o (u + d') k 


sin(7r<5) r(fc)r(<J-fc + l) 




sin(7T(5) -f—j' \k J r(fe)r(<5 — k + 1) 

sm(7n5) ^ \kJ \k — 1 J 

For the -1-st moment, we obtain 

f_! = -^R 2 r(i + <y)r(i - 8)e s P ( 1 - P ) 5 ~\ P <i, 


M_i = exp(C l # fl p(l — P ) S *) 


= M 


-(i -pY 

i 


, p < 1. 


and thus 
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